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t Abstract 

It is shown that for the spherically-symmetric and static systems the hypotheses 
„Cj , posed by Yang and Radinschi and by Vagenas can be related to the particular distri- 

bution of the source. Simple proofs are given and a number of examples are discussed 
with the special emphasis put on the quantum corrected Schwarzschild black hole. 
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Recent interest in the pseudotensors describing the energy-momentum of gravity has 
$— i | been stimulated to a large extent by the understanding that the freedom in defining the 

pseudotensors is associated with the freedom in the choice of the Hamiltonian boundary 
terms pQ. Thus, each pseudotensor reflects particular physical situation, and, being in fact 
quasilocal [2J, is acceptable. This removed (at least partially) the anathema pronounced 
by a number of authors (see for example Ref. 3) and caused that the energy momentum 
pseudotensors returned from the exile and have settled themselves in the main stream of 
the theory. There are quite a number of various constructions, each of noble parentage: 
Einstein0], Landau-Lifshitz [5], Weinberg [6], M0ller[7], Bergman [8] and Papapetrou [9], to 
name a few. Each has been a subject of detailed studies and a source of interesting and 
important results. The subjects covered range form the black hole physics to gravitational 
waves. (See, for example, Refs. 10-18 and references cited therein). Despite extensive work 
■ done so far there is still room for new findings and fresh insights. 

Some time ago Radinschi and Yang [17] observed that for the Reissner-Nordstrom and 
the regular black holes [HI (201 [23 [22] the difference between the energy calculated with the 
aid of the Einstein pseudotensor, Ee , and this obtained employing the M0ller prescription, 
Em , can be related to the energy density, /?, of the matter fields 

E E -E M ~ r 3 p. (1) 

Later Vagenas [H] found that for the Reissner-Nordstrom (RN) and the RN dyadosphere 
[23] there is an interesting relation between the coefficients of the expansion of Ee and Em 
in the inverse powers of r. This has been partially confirmed by the behaviour of a few first 
coefficients of the expansion of the regular black hole solution. He hypothesized that if 
and Pk are, respectively, the expansion coefficients of Ee and Em, then 

a k = j^P, (2) 
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It is the aim of this letter to demonstrate how these observations and hypotheses can be 
proven and understood in relation with the source term of the Einstein field equations, (i. c., 
the stress-energy tensor of the matter fields) and to analyse their limitations. The discussion 
will be illustrated with a few examples. As our demonstration heavily relics on the particular 
representation of the line element describing static and spherically-symmetric configuration 
we shall start with a few basic facts. 

Let us consider the general spherically-symmetric, static and asymptotically flat geometry 
described by the line element (our conventions follow these of MTW) 

ds 2 = _ e 2V>M A _ ^ d e + (l - 1 dr 2 + r 2 (d9 2 + sin 2 9d0 2 ) , (3) 

where M (r) and ip (r) are two functions that are to be determined from the Einstein field 
equations, supplemented by the physically motivated boundary conditions. This represen- 
tation of the line element has proven to be useful in various contexts as the Einstein field 
equations have particularly transparent form [24] . Moreover, let us assume that the (total) 
covariantly conserved stress-energy tensor of the material fields can be written in the form 
dictated by the symmetries of the problem 

T£ = diag [a (r) ,b(r),c(r),c (r)] . (4) 

The (q) component of the Einstein equations reads 

2 dM (r) 



r 2 dr 

whereas their (I) component can be written as 



8™ (r) , (5) 



2 dM(r) 2 / 2M(r)\ , 

W - 1 — V (r) = Snb (r) . (6) 



r 2 dr r 

In the latter we shall employ the boundary conditions in the form 

M(oo)=M, VM = 0, (7) 

where M. is the total mass of the system as seen by a distant observer. On the other hand, 
for the black hole configurations, one can also take 

M{r + ) = r -±, V(oo)=0, (8) 

where r + is the location of the event horizon, defined here as the outermost trapped surface. 
From (O one has 

dM(r)_ „ , 



dr 

whereas from ([5]) and ^ one obtains 



4nr 2 a (r) , (9) 



^1 ^ inr ^Ll^M ( 10 ) 

dr 1 2M(r) ■ y 1V > 

r 

Now, let us return to the general form of the line element and calculate the Einstein 
pseudotensor The Einstein pseudotensor, 0^, can easily be constructed in the Cartesian- 
type coordinates by differentiation of the Freud superpotential, V^ v , 
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where 
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Vr = —9xr— [9 (iTsT - sTg™)] ■ (12) 
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The 0^ does not exhaust the list of interesting and important pseudotensors. Indeed, 
there are half a dozen other objects (sometime referred to as complexes) that are useful in 
various circumstances. One of the most studied, however, is the pseudotensor constructed by 
M0ller in the late fifties of the last century. One of the advantages of the M0ller pseudotensor, 
7jf*, is that it can be calculated in any coordinate system. Its construction proceeds in two 
steps: first, one has to calculate the superpotential 

^ A = -V=a (-Lg** ~ tL*>*) (Th*, (is) 



v dx K dx 
and, subsequently, the energy-momentum pseudotensor 
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Due to the antisymmetry of the superpotentials the pseudotensors (fTTj) and (|14j) are 
divergence-free. The energy and momentum are given by 

P t = JJJ T°dx 1 dx 2 dx 3 (15) 

and 

E = III rgdx 1 dx 2 dx 3 , (16) 



respectively, where r/f stands either for O^ or T^. In the latter we shall differentiate between 
them simply by inserting subscripts E or M. 

The energy calculated with the aid of the M0ller energy-momentum pseudotensor for the 
geometry (j3]) assumes the simple form 

E M = [M (r) - rM' (r) + r 2 ^' (r) - 2rM (r) tji' (r)] . (17) 

On the other hand, the calculations of the Einstein pseudotensor is slightly more involved. In- 
deed, as have been said earlier one has to introduce the standard quasi-Cartesian coordinates 
{x, y, z} and rewrite the line element ^ in the form 

ds 2 = „ e 2^(r) ^ _ ^ + (dx 2 + dy 2 + dz 2 ) 

2M(r) 



(xdx + ydy + zdz) 2 . (18) 



From the definition of the Einstein pseudotensor, after some algebra, one obtains remarkably 
simple result 

E E =e^ r) M{r). (19) 

The calculations have been carried out with the aid of Maxima package, and in the advent 
of the computer algebra systems, this problem should be regarded as a five-finger exercise. 
Identical results can easily be obtained from the formulas constructed by Virbhadra 
and Xulu for the Einstein and M0ller pseudotensors, respectively, in the most general 
nonstatic and spherically-symmetric geometry. 

Now, let us assume that a (r) falls sufficiently rapidly as r — > oo to make the integral (J9j) 
finite and a(r) = b (r) . Making use of the boundary conditions ([7]) one has ip (r) = and 

M (r) = M + 4tt / r 2 a (r) dr = M + M (r) . (20) 
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On the other hand, the boundary conditions ((5J) yield 



M (r) = — + 4tt / r 2 a (r) = — + M (r) . (21) 
2 J r+ 2 

For the simplified configuration the M0ller energy reads 

E M = M (r) - rM' (r) , (22) 

whereas the Einstein energy coincides with the function M (r) as can easily be seen from 
Eq.®: 

E E = M (r) . (23) 

Now, we are ready to demonstrate a simple relation existing between the expansions (in 
terms of r" 1 ) of the Einstein and M0ller energy. Following Vagenas, let us assume that 

oo 

M (r) = ^a fe r- fe , (24) 
fc=i 

where a k are numerical coefficients. It should be noted that there are quite a number of 
solutions that belong to that class. Indeed, the Reissner-Nordstrom, the regular black holes 
or the geometry of the dyadosphere of RN black hole may serve as the simple examples of 
particular realizations of (|24|) with 

O 2 

<*l = ~, a k = for k > 2, (25) 



and 

Q 2 Q 4 
at = fl 2 = a 3 = a i — 0) a 5 = a k — for fc > 6, (27) 

respectively. We do not attempt, of course, to list all the solutions of this type here. 
It can easily be demonstrated that if M (r) is given by ([24| then 

oo 

E E =M + Y, a kr- k (28) 

k—l 

and 

oo oo 

E M =M + ^2(k+l) a k r- k = M + £ (3 k r- k , (29) 
fc=i fc=i 

and, consequently, 

^ = ^A. (30) 

The above equation is precisely the relations found by Vagenas |18] . He correctly hypothesized 
that it should have a wider domain of applicability than two or three special cases considered 
in Ref. 18. It should be emphasized, however, that in this demonstration the crucial role is 
played by the special form of the stress-energy tensor. 

On the other hand, still working with the simplified stress-energy tensors of the matter 
fields one has 

AE = E M -E e = -rM' (r) 

= 47rr 3 a(r) = 47rr 3 T °, (31) 
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where we have used Eq. @. This result for particular configurations has been obtained by 
Yang and Radinschi in Ref. 17. Finally observe, that for a(r) ^ b (r) the relations (|30[) and 
(f3"Tj) do not hold in general. 

Equipped with the results (|22|) and (|23|) let us consider a slightly more complicated 
example of the quantum corrected black hole being a solution of the semiclassical Einstein 
field equations. The idea that lies behind the calculations is simple. One has to construct the 
stress-energy tensor for a general background ^ and subsequently solve (self-consistently) 
the differential equations with the boundary conditions |(7J or ©. 

The stress-energy tensor of the quantized massive scalar, spinor and vector fields can 
be obtained from the one-loop renormalized effective action, Wr, in a standard way, i.e., 
by a functional differentiation of Wr with respect to the metric tensor. It gives rise to the 
quantum-corrected black hole described by the line element The most general expression 
describing the approximate renormalized stress-energy tensor of the quantized massive scalar, 
spinor and vector fields have been constructed in Refs. 27 and 28. These formulas generalize 
the earlier results obtained by Frolov and Zel'nikov [29l [3Ql [31] (valid for the Ricci-flat 
geometries) and those of Ref. 32, where the stress-energy tensor has been calculated for the 
massive scalar fields in the general static, spherically-symmetric geometries. 

In what follows we shall employ the boundary conditions ([7]), i.e., we shall parametrize 
the solution by the total mass of the system (consisitng of both classical and quantum parts) 
rather than the bare black hole mass. The resulting differential equations are, in general, too 
complicated to be solved exactly and consequently one has to refer to approximations. For- 
tunately, one can easily adopt an approach in which the stress-energy tensor of the quantized 
fields can be regarded as a small perturbation. Indeed, this can be done simply because the 
total stress-energy tensor of the matter fields can be written as 

T£ = diag [a (r) + £ (T*) « , b (r) + s(T^ n ,c (r) + e <I*>« , c (r) + e{T*)&] , (32) 

where (T^)i s J n is the approximate renormalized stress-energy tensor of the quantized massive 
scalar (s = 0), spinor (s = 1/2) and vector (s = 1) fields, and e is the auxiliary parameter 
that helps to keep track of various terms in complicated expansions. (One should put e = 1 
in a final stage of calculations). 

To simplify discussion let us assume that for the classical sources the condition a (r) = b (r) 
holds. This allows us to represent the functions M (r) and ip (r) , with a small abuse of 
notation, in the forrrj^ 

M (r) = M (r) + eM x if) + O (e 2 ) and ip (r) = eipi (r) + O (e 2 ) . (33) 

Now, restricting to the first-order perturbations one has the following system of the differen- 
tial equations for Mi (r) and ipi (r): 

^ = -4- 2 (T*)(i (34) 

and 

- 47rr 2M (r) • 

r 

The energy calculated from the Einstein pseudotensor is 

E E = M (r) + [M (r) Vi (r) + M[ (f)\ e + O (e 2 ) (36) 



1 The function Mq (r) in Eg. 1331 must not be confused with the one defined in Eq. 1201 
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whereas the M0ller energy is given by 



E M = M (r) - rM' Q (r) + [M (r) fa (r) + Mf (r) - r^i (r) M' Q 

-rM[ (r) + rVi (r) - 2rM (r) $ (r)] s + O (e 2 ) ; (37) 

and the difference AE is 

AE = -rM' Q (r) - [fa (r) Mq + rAf{ (r) - r 2 i/4 (r) + 2rAf (r) (r)] e + (e 2 ) . (38) 

In the simplest case, the classical sources are absent and the zeroth-order solution reduces 
to the Schwarzschild line element with M (r) — M. (Note that M (r) = M and Mi (r) = 
as r — > oo). Although the renormalized stress-energy tensor in the Schwarzschild spacetime 
has been calculated by Frolov and Zel'nikov [55] in the present calculations we have used 
general formulas of Refs. 27 and 28. After some algebra, for the quantized massive scalar 
field (s = 0) in a large mass limit with arbitrary curvature coupling, one has 



M 2 

M^(r)=M 



19 313 
(lLM-6r)jj+— r-— M 
56 504 



157rm 2 r 6 
and 



(39) 



( r ) = _„ _ _ . (40) 

Here m is a mass of the field whereas rj = £ — 1/6 is the numerical parameter and its two 
particularly appealing values are £ = for the minimal and £ = 1/6 for the conformal 
coupling. Similarly, for the quantized massive spinor the first-order solution is simply given 

by 

M^nr)=M--J^Ur- l 4M], (41) 



140vrm 2 r 6 \ 27 
<A(V2) (r) - (42) 



whereas for the quantized massive vector field one has 



(r)=M + 55^53 ( 37r - ) , (43) 



2807rm 2 r 6 V 9 

^ (r) = 131M *. (44) 
v K ' 8407rm 2 r 6 v ; 

Eqs. ((36| . ([37]) and (|39M4II are, of course, sufficient to calculate the Einstein and M0ller 
energy. On the other hand, combining (j3"5|) , (|34[) and (|35p . for a problem on hand, for AE 
one has 

A J B = 4^r 3 (r;)(i. (45) 

It is an interesting result, which clearly shows that the difference between the M0ller and 
Einstein energy is related to the radial pressure of the quantized fields. It can be contrasted 
with the Radinschi-Yang formula. One can ascribe this behaviour to ip (r) ^ 0, which, in 
turn, is a consequence of the form of the stress-energy tensor. On the other hand, for the 
class of the renormalized stress-energy tensors of the quantized fields in a large mass limit in 
the Bertotti- Robinson (or more generally AdS2 x S%) geometries one has ip (r) — 0. It should 
be noted that the formula ([45|) may be used as the useful check of the calculations. 
Now, making use of either (|3"8|) or (jl5|) . one readily obtains 



AE^= M 



2 r 6 



47rr 3 (T;>(° e ) n , (46) 
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ajbMO = ( I- - I.m ) = W(t;)(v 2 ) (47 ) 



^ 2 r 6 \ 20 90 

and 



A£ (1) = ( 1 4-^M)= . (48 ) 



^ 40 12 

It should be emphasized once more, that (T")rJn calculated for the zeroth-order solution 
of the semiclassical Einstein held equations coincides with the stress energy tensor of the 
massive fields constructed in the Schwarzschild spacetime parametrized by the mass M.. 

To this end, let us calculate the Einstein energy as the M0ller energy can easily be obtained 
from the definition of AE^ and (|46M8p . Inserting the functions (r) and ip^ ( r ) as 
given by (|39M4p into (|36p (or, equivalently, into (fT9|) and retaining the terms which are linear 
in e) one obtains 



JR 157rm 2 r 6 



59 19 

—M--r+(6r-AM)v 



(49) 



and 



41/2) =M + ( _L M _ J_ r 1 (50) 

E irm 2 r e \ 378 140 ' 



p (i) . . M /218 



E *- M -m^[ir M - 37r i- { ^ 



Finally let us return to the boundary conditions (© and observe that in this case the 
solution is parametrized by the exact location of the event horizon, r+. The total mass of 
the system is, therefore, given by 

M = lim (M 0(r) + eM x (r)) , (52) 

r — >oc 

and it differs from the horizon defined mass, Mh . As the calculations in this parametrization 
closely follow the general pattern described above, we shall not display them here. 

Our next example is to certain extent related to the previous one: we shall consider 
the pure gravity action supplemented with the higher curvature terms and analyse the 
Schwarzschild solutions influenced by all time-reversal-invariant operators of dimension six 
[3"3"j . In doing so we shall assume that the coefficients of the higher curvature part of the 
action are arbitrary rather than inspired by a particular theory. The action functional has, 
therefore, the form 

S = J d A xg l/2 R + S l + S 2 i (53) 

where S\ is the action functional of the quadratic gravity (that does not influence the prob- 
lem) and S 2 is given by 



S 2 = J d 4 xg 1/2 (cxR^R^R^^ + c 2 R 
+czR l _ lv R^ a ^ 1 R v a/ 3 7 + CiRR l _ lv \ K R^ x ' i 



+ c 5 R^R XK R^ K + c 6 R^R^ x R x v + c 7 R 3 + c 8 RR^R^) . (54) 
For spherically-symmetric and static systems the perturbative solution is given by ([3]) with 

967T 

ip{r) = 7-M 2 e (18ci + 6c 2 + 71c 3 + 28c 4 ) (55) 
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and 



M (r) = M g e (12ci + 3c 2 + 4c 3 + 16c 4 ) 



+ 32 ^ 3 £ (98ci + 25c 2 + 33c 3 + 132c 4 ) . (56) 

Elementary calculation shows that the difference between Einstein and M0ller complexes is 
given by 

AE = -rM[ (r) + r 2 ^ (r) (l - . (57) 

Making use of (f5H|) and ([ST]) , one arrives at 

1447T.M 2 

AE = - (16ci + 9c 2 + 8c 3 + 32c 4 ) 

- 192 ^ 3 (10ci + llc 2 + 9c 3 + 36c 4 ) . (58) 

It can easily be demonstrated by an independent calculation that 

S r r = AE/8nr 3 , (59) 

where 

S ab =4j- 2 ^. (60) 
9 1/2 og ab 

Treating S ab as the effective stress-energy tensor, one concludes that AE is proportional to 
the radial pressure, as expected. 

The results obtained in this paper clearly shows that it is necessary to explain why 
various pseudotensors are related in a specific way. We suspect that there are other, quite 
simple but interesting relations which can easily be explained with the aid of the Einstein 
equations and careful examination of the role played by their source terms. Moreover, it 
would be interesting to consider other pseudotensors for the whole classes of geometries and 
relate them with the specific choices of the boundary terms. Finally we observe that the 
evident direction for further research is generalization of the semi-classical calculations to 
more complicated settings, such as the Reissner-Nordstrom spacetime or the geometry of 
the regular black holes. These issues are under active consideration and the results will be 
published elsewhere. 
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